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Abstract. In this paper we report a study (see Mastrobuono-Battisti et al., 2012) about the 
formation and characteristics of the tidal tails around Palomar 5 along its orbit in the Milky 
Way potential, by means of direct N-body simulations and simplified numerical models. 
Unlike previous findings, we are able to reproduce the substructures observed in the stellar 
streams of this cluster, without including any lumpiness in the dark matter halo. We show 
that overdensities similar to those observed in Palomar 5 can be reproduced by the epicyclic 
motion of stars along its tails, i.e. a simple local accumulation of orbits of stars that escaped 
from the cluster with very similar positions and velocities. This process is able to form 
stellar clumps at distances of several kiloparsecs from the cluster, so it is not a phenomenon 
confined to the inner part of Palomar 5's tails, as previously suggested. 
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1. Introduction 

Palomar 5 (Pal 5) is an outer Milky Way glob- 
ular cluster (GC) which presents some pecu- 
liar characteristics: a very low mass, M tol = 
5 xlO 3 M , a large core radius, r c = 24.3 pc, 
and low central concentration, c=0.66 (see 
Odenkirchen et al., 2002). Moreover, it is sur- 
rounded by two massive tidal tails emanat- 
ing from opposite sides of the cluster and 
populated by stars escaped from the system 
(Odenkirchen et al., 2001, 2003; Grillmair & 
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Dionatos , 2006; Jordi & Grebel, 2010). With 
an overall detected extension of 22° on the sky, 
corresponding to a projected spatial length of 
more than 10 kpc, Pal 5's tails are so elon- 
gated and massive that they contain more stel- 
lar mass than the one currently estimated to be 
in the system itself (Odenkirchen et al., 2003). 
Pal 5's streams are characterized by the pres- 
ence of inhomogeneities: stellar density gaps 
(underdense regions) and clumps (overdense 
regions) are particularly visible in the trail- 
ing stream, with the most evident and massive 
overdensities found between 100 and 120 ar- 
cmin from the cluster center (Odenkirchen et 
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al., 2003). As Capuzzo-Dolcetta et al. (2005) 
and Di Matteo al. (2005) showed, the presence 
of similar substructures in globular cluster tidal 
tails can be related to kinematic effects and lo- 
cal decelerations in the motion of stars, later 
identified by Kiipper et al. (2008, 2010, 2012) 
and Lane et al. (2012) as "epicyclic cusps", 

1. e. regions where stars escaping from the clus- 
ter slow down in their epicyclic motion, as 
seen in the cluster reference frame. However, 
while all these works were successful in show- 
ing that complex morphological substructures 
can be formed in the tidal streams of globu- 
lar clusters, until now no simulation or model 
has been able to reproduce convincingly the 
only stellar clumps robustly observed in the 
tails of a Galactic GC: those of Pal 5. Even 
in the most complete numerical study of Pal 
5's stellar streams realized in the last years, 
Dehnen al. (2004) have not succeeded in re- 
producing the clumpy nature of the stream. 
This lead them to propose that the observed 
clumps may be the effect of Galactic substruc- 
tures (as giant molecular clouds or dark matter 
subhalos) not accounted to in their simulations. 
The hypothesis that the gravitational effects of 
dark subhalos can generate a complex mor- 
phology in tidal tails, with overdensities and 
gaps, has been recently confirmed by Yoon et 
al. (201 1). However, in the results presented by 
Mastrobuono-Battisti et al. (2012) and summa- 
rized in this paper, we show that the inhomo- 
geneous tails of Pal 5 can be reproduced also 
without invoking any lumpy dark halo, just as 
an effect of epicyclic motion of stars. 

2. Models and initial conditions 

2.1. N-body simulations and restricted 
three-body methods 

We used N-body, direct summation simula- 
tions to investigate the formation and charac- 
teristics of Pal 5's tidal tails (see Mastrobuono- 
Battisti et al., 2012, for more details). To run 
the simulations we used a modified version of 
NBSymple (Capuzzo-Dolcetta et al., 2011), a 
high performance direct A^-body code imple- 
mented on a hybrid CPU+GPU platform. In 
this code, the effect of the external galactic 



field is taken into account by mean of an ana- 
lytical representation of its gravitational poten- 
tial (see Sect. 2.2). To smooth the mutual in- 
teraction between particles we adopted a soft- 
ening length, s = 0.03 pc. The time integra- 
tion of the particles trajectories was done using 
the second-order leapfrog method. The average 
relative error per time step was less than 10~ 13 
for the cluster both isolated or in the Galactic 
potential. The result of the N-body simulations 
were compared to simplified, restricted three- 
body models, where we used the same initial 
conditions adopted to run the simulation but, at 
each time step, instead of considering the mu- 
tual interactions, we evaluated the interaction 
with the global gravitational potential of the 
cluster and with the Galactic potential for each 
particle. Finally, similar to Lane et al. (2012), 
we produced streaklines to visualize in a more 
clear way the spatial distribution, at the present 
time, of stars that have escaped from Pal 5 in 
the last Gyrs. To do so, we integrated the orbit 
of a point mass representing Pal 5 barycenter, 
and, at each time step, released two particles at 
distance r max = 120 pc from it, along the in- 
stantaneous Galactic center and anticenter di- 
rections, and with a velocity equal to that of 
the cluster. The motion of these escapers is in- 
tegrated in time, taking both the effects of the 
GC and Galactic potential into account. 

2.2. Globular cluster and Galaxy models 

The Galactic potential is from Allen & 
Santillan (1991), whose model consists of a 
spherical central bulge and a flattened disk, 
plus a massive smooth spherical halo. For the 
GC we adopted the parameters of the "model 
A" described by Dehnen al. (2004) which con- 
sists of a single-mass King model (King, 1966) 
with tidal radius r, = 56 pc, Wo = 2.75, and 
M = 2xlO 4 M ; we used = 15360 particles 
for the N-body representation of this model. 

To obtain the orbital initial conditions, we 
integrated a test particle with Pal 5's present 
position and velocity backward in time for the 
same time interval as chosen by Dehnen al. 
(2004), i.e. 2.95 Gyr, in the Galactic potential. 
Then we translated the positions and velocities 
of the King model to coincide with those of 
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Fig. 1. Upper panel: The trailing tail at the end 
of the N-body simulation. The contour lines in 
the tails refer to 740, 2200, 4800, 6300 and 
7400 M Q /deg 2 . Bottom panel: Linear density 
of the trailing tail shown in the upper panel, as 
a function of the right ascension a. The clus- 
ter is on the right of the panel. Error bars are 
Poisson errors. From Mastrobuono-Battisti et 
al. (2012). 

the cluster barycenter 2.95 Gyr ago, and we in- 
tegrated the whole cluster 2.95 Gyr forward in 
time. 

3. Results and discussion 

During its evolution, the simulated cluster has 
lost the most of its initial mass (88%), now 
redistributed into two long and narrow tidal 
tails. This tails are only apparently smooth: as 
it becomes evident plotting the cluster isoden- 
sity contours, their structure is actually com- 
plex and the stellar distribution is clearly inho- 
mogeneous. In particular, as shown in the up- 
per panel of Fig. 1, the densest substructures 
are all localized in the trailing tail. In this tail, 
our N-body simulation predicts there are two 
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Fig. 2. The trailing tail of the cluster in the re- 
stricted three-body model. The contour lines 
in the tails refer to 470, 940, 1410, 1970, 
2500 and 3300 M /deg 2 . From Mastrobuono- 
Battisti et al. (2012). 



prominent density enhancements: the first lo- 
cated very close to the cluster (a < 230°) and 
the second starting at ff=230.5° and ending at 
a=232°. The second region shows an under- 
dense region in it, whose extension is ~ 0.5°. 
All these characteristics, at the same location, 
are also found in the observed streams (see 
Fig. 3 in Odenkirchen et al., 2003). Not only 
does the position of these clumps closely re- 
semble the location of those observed in Pal 
5's tails, but their linear densities (bottompanel 
of Fig. 1) - 2-3 times above the density of 
the surrounding streams - and the presence 
of underdense regions, with sizes of 0.5°-l°, 
also agree with those measured for Pal 5 (see 
Fig. 4, Odenkirchen et al., 2003). As shown 
in Fig. 2, these substructures are also visible 
in the isodensity contours of the simplified nu- 
merical model, at few degrees from the cluster 
center. These regions of over- and underden- 
sities are due to the epicyclic motion of stars 
escaping from the cluster: stars lost at differ- 
ent times redistribute along the tail following 
a complex path, as shown by the streaklines in 
Fig. 3. This kinematic process, described and 
studied in a number of papers (Kiipper et al., 
2008, 2010, 2012; Lane et al., 2012), is ap- 
plied here for the first time to reproduce ob- 
served stellar streams, including both the po- 
sition and the intensity of the observed stellar 
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Fig. 3. Streaklines which show the positions, 
at the present time, of stars that have escaped 
from Pal 5 in the last 2.95 Gyrs. Different col- 
ors correspond to stars lost during different 
time intervals of 0.5 Gyr. From Mastrobuono- 
Battisti et al. (2012). 



inhomogeneities. The epicyclic motion is thus 
able to form over- and underdensities in sev- 
eral regions of the tails, at a variety of distances 
from Pal 5. If the orbit is integrated over longer 
times (~8 Gyr), streaklines indeed show that 
epicycle loops also form at distances of several 
kpc from the center. Moreover, since these sub- 
structures are also reproduced in models where 
the mutual gravity is neglected, it is possible to 
rule out Jeans instability as a possible mech- 
anism responsible for their formation, as sug- 
gested by Quillen & Comparetta (2010). 



4. Conclusions 

We have studied the formation and character- 
istics of the tidal tails of the GC Pal 5 along 
its orbit in a smooth Milky Way potential by 
means of N-body simulations and simplified 
numerical models. For the first time, we were 
able to reproduce the observed inhomogeneous 
structure of Pal5's tidal tails using N-body sim- 
ulations, without including any lumpiness in 
the Galactic halo. The epicyclic motion of stars 
along the tails is the main mechanism respon- 
sible for the formation of clumps. This result 
must be taken into account when using streams 
inhomogeneities to evaluate the granularity of 
the Milky Way dark halo. 



Acknowledgments 

This work was carried out under the HPC- 
EUROPA2 project (project number: 228398) 
with the support of the European Commission 
Capacities Area-Research Infrastructures 
Initiative. 

References 

Allen, C. & Santillan, A. 1991, RMxAA, 22, 
255 

Capuzzo-Dolcetta, R., Di Matteo, P. & 

Miocchi, P. 2005, AJ, 129, 1906 
Capuzzo-Dolcetta, R., Mastrobuono-Battisti, 

A. & Maschietti, D. 201 1, NewA, 16, 284 
Dehnen, W., Odenkirchen, M., Grebel, E. K. & 

Rix, H.-W. 2004, AK, 127, 2753 
Di Matteo, P., Capuzzo Dolcetta, R. & 

Miocchi, P. 2005, Celestial Mechanics and 

Dynamical Astronomy, 91, 59 
Grillmair, C. J. & Dionatos, O. 2006, ApJL, 

641,37 

Jordi, K. & Grebel, E. K. 2010, A&A, 522, 71 

King, I. R. 1966, AJ, 71,64 

Kupper, A. H. W., MacLeod, A. & Heggie, D. 

C. 2008, MNRAS, 387, 1248 
Kupper, A. H. W., Kroupa, P., Baumgardt, H. 

& Heggie, D. C. 2010, MNRAS, 401, 105 
Kupper, A. H. W., Lane, R. R., & Heggie, D. C. 

2012, MNRAS, 420, 2007 
Lane, R. R., Kupper, A. H. W. & Heggie, D. C. 

2012, MNRAS 423, 2845 
Mastrobuono-Battisti, A. , Di Matteo, P., 

Montuori, M., Haywood, M. 2012, A&A, 

547, L7 

Odenkirchen, M., Grebel, E. K, Rockosi, 
C. M., et al. 2001, ApJL, 548, 165 

Odenkirchen, M., Grebel, E. K., Dehnen, 
W., Rix, H.-W., Cudworth K. M. 2002, AJ, 
124, 1497 

Odenkirchen, M., Grebel, E. K, Dehnen, W., 
Rix, H.-W., Yanny, B. et al. 2003, AJ, 126, 
2385 

Quillen, A. C. & Comparetta 2010, MNRAS 

submitted; astro-ph/1002.4870 
Yoon, J. H, Johnston, K. V., Hogg, D. W. 201 1, 

ApJ,731,58 



